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We present methods for the direct characterization of quantum dynamics (DCQD) in which both the princi¬ 
pal and ancilla systems undergo noisy processes. Using a concatenated error detection code, we discriminate 
between located and unlocated errors on the principal system in what amounts to filtering of ancilla noise. The 
example of composite noise involving amplitude damping and depolarizing channels is used to demonstrate the 
method, while we find the rate of noise filtering is more generally dependent on code distance. Our results indi¬ 
cate the accuracy of quantum process characterization can be greatly improved while remaining within reach of 
current experimental capabilities. 


I. INTRODUCTION 

Decoherence and noisy dynamics in open quantum systems 
are major hurdles to the realization of working quantum com¬ 
puters and practical quantum communication devices [OQ • In¬ 
formation encoded in a controlled quantum system will leak 
into its surrounding environment when there is coupling be¬ 
tween the two systems, resulting in shorter qubit lifetimes 
and lower gate fidelities. It is necessary to characterize these 
quantum dynamics in order to better understand the sources 
of system-environment coupling. Characterization of the dy¬ 
namical process can then be used to mitigate sources of noise 
and improve qubit coherence times. 

Given the density matrix p for a d-dimensional system, the 
completely-positive quantum process 

d 2 -1 

£{p) = 'y ] EqpE^ = y ] XmnFmpFn (I) 

a m,n =0 

may be described in terms of its Kraus operators {E a } or the 
Hermitian process matrix x defined with respect to the opera¬ 
tor basis {F m }. Experimental characterization of the process 
matrix x provides a concrete representation of £ that can be 
used to study and refine system behavior. In standard quan¬ 
tum process tomography (SQPT), measurements characteriz¬ 
ing the state £(p) are used to reconstruct the process matrix 
by inverting Eq. 0 over a complete set of input states m- 
[61. Ancilla-assisted process tomography (AAPT) performs a 
similar inversion using fewer input states by exploiting corre¬ 
lations between the principal system (P) and an ancilla sys¬ 
tem (A) isolated from non-trivial quantum process Q. Such 
a composite process can be written as £p ® 1 a(p) where the 
subscripts indicate which processes occur on each subsystem. 

In contrast, direct-characterization of quantum dynamics 
(DCQD) avoids inverting Eq. 0 by measuring the process 
elements Xmn directly Glia. DCQD techniques have re¬ 
cently been applied to characterize trapped ion (H and hyper¬ 
entangled photon CD dynamics. Like AAPT, the princi¬ 
pal and ancilla subsystems are initially entangled in a probe 
state before being subjected to non-trivial and trivial quan¬ 
tum processes respectively. Interestingly, DCQD probe states 
can be described as the codewords of a quantum error cor¬ 
rection (QEC) code. In this framework, a quantum process 
maps the joint system probe state either within or outside 


the codespace. Processes mapping the probe state outside 
the original codespace are detected and characterized by their 
error syndrome, i.e., the measured eigenvalues of each QEC 
code generators. Syndrome frequencies derived from an en¬ 
semble of stabilizer measurements are sufficient to directly 
characterize the underlying process matrix x- 

The DCQD framework shows that the mathematical tools 
developed for QEC can be leveraged for process character¬ 
ization. In particular, code design plays an important role 
in probing the process matrix EH. Recent extensions to 
DCQD involve generalized characterization codes which also 
encode logical quantum information [12, 131. This offers the 
ability to characterize processes occurring during an arbitrary 
quantum computation. Complimentary works, from a QEC 
perspective, have shown that syndrome data generated by er¬ 
ror correction protocols can be used for noisy parameter esti¬ 
mation OHS, with recent experiments involving stabilizer 
QEC circuits over 9 and 4 qubits are prime candidates for 
these types of characterization methods [ 17, 18|. 

Despite advances within the DCQD paradigm, a significant 
and persistent drawback in all existing schemes is the require¬ 
ment that the ancilla system be perfectly noiseless. This as¬ 
sumption is necessary for correctly interpreting the measured 
syndromes in the context of Eq. 0. Noisy ancilla lead to spu¬ 
rious data that corrupts the process tomography and adds er¬ 
rors to the process matrix. However, noise is certainly present 
in any realistic experiment and it is important to ask how 
QEC-based process characterization can be extended to in¬ 
clude noisy ancilla. 

We address the use of noisy ancilla for process characteriza¬ 
tion by introducing a new class of quantum process codes that 
remove the requirement of noise-free ancilla. Our approach 
is based on concatenated encoding of the ancilla system us¬ 
ing a second quantum error detection code. We show that 
by monitoring syndrome values of the composite code mea¬ 
surements of the principal system that have been corrupted 
by ancilla noise can be filtered out. By removing measure¬ 
ments attributed to noisy ancilla, we generate a higher fidelity 
construction of the process matrix than possible with direct 
characterization alone. We also examine the question of effi¬ 
ciency, which we define as a tradeoff between the syndromes 
collected and the accuracy of the process characterization. 

The remainder of the paper is organized as follows: in 
Sec. [II] we review notation for stabilizer QEC codes and out- 
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line the conventional DCQD procedure for constructing the 
process matrix before introducing a concatenated six-qubit 
code used to characterized the dynamics of a two-qubit princi¬ 
pal system in the presence of full system noise. This includes 
a discussion of how ancilla error detection is used to filter to¬ 
mographic information prior to characterization. In Sec. [Ill] 
we present a numerical case study of an amplitude damp¬ 
ing channel on various codes with and without depolarizing 
noise affecting the ancilla subsystem. We discuss results of 
our simulation, the degree to which the code faithfully charac¬ 
terizes dynamics on the principal system, and the probability 
that high weight errors, which can pass through our concate¬ 
nated error filter thus corrupting the tomographic data, occur 
in Sec. IV Our conclusions and discussion appear in Sec.fV] 


II. CHARACTERIZATION WITH NOISY ANCILLA 

An [[n, d]\ stabilizer code maps k logical qubits onto n 

qubits with a distance d between distinct codewords EES). Let 
S = {gir ’ ' 1 9r) denote an Abelian stabilizer group whose 
r = n — k generators are drawn from the n-qubit Pauli group, 
i.e., gi G V n . These stabilizer generators define a set of com¬ 
muting observables called the syndrome that partitions the 
n-qubit Hilbert space into a set of mutually orthogonal sub¬ 
spaces, each encoding k qubits. The i-th subspace corre¬ 
sponds to a syndrome eigenvalue e^. We will represent each 
syndrome as a string of classical bits such that = 0 or 
1, respectively, for the +1 or —1 eigenstate of the generator 
gj. In this notation, the logical codespace J #o corresponds 
to the trivial syndrome eo generated by the stabilizer group 
S = {seP n :s\il>) = \il>), V|^)G^o}. 

Errors due to a set of operators E are said to be correctable 
if the QEC condition 

(i\EtE b \j) = CabSij (2) 

is satisfied for all E a ,Ei) G E, where |i), \j) are orthonor¬ 
mal basis vectors spanning J #o and C a b is a Hermitian ma¬ 
trix. In particular, a correctable error E a maps states in the 
codespace ^ to another codespace. Upon measuring a syn¬ 
drome value ei , a state |^) is projected into the subspace 
For purposes of error correction, the syndrome dictates what 
recovery operation should be applied to return the state to the 
logical codespace 

In the context of quantum process characterization the op¬ 
erator set E represents a basis for the dynamical processes by 
which the encoded state evolves. Instead of correcting E, the 
goal of process characterization is to unambiguously detect 
these operations. As such, there is a significant difference be¬ 
tween how QEC and process characterization are affected by 
undetectable errors. In particular, operators commuting with 
the stabilizer group belong to the normalizer group A f(S), 
where elements of the normalizer act as logical operators 
that map one codestate to another or stabilize the state (since 
S G A Elements of A f(S) yield trivial syndromes and 
lead to logical errors because no recovery operation is applied. 
For process characterization, however, undetectable stabilizer 
and normalizer operators lead to faulty tomographic data as 
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FIG. 1. Outline of a quantum circuit used to directly characterize the 
dynamics of a principal system P. The principal system is first entan¬ 
gled with the ancilla system A before being subject to some noisy 
dynamics. The process matrix is then constructed from an ensem¬ 
ble of stabilizer measurements (for details see Sec. |II A| and Fig. [2j. 
Panel (a) illustrates a conventional DCQD circuit, where the P and A 
are maximally entangled in a Bell state |<I> + ) and characterization as¬ 
sumes a noiseless ancilla. Panel (b) illustrates a concatenated DCQD 
circuit, with an initial entangled state |0) (Eq. [8j, which supports 
errors on the ancilla subsystem. 


they cannot be distinguished from the identity operation. As 
seen below, concatenated ancilla codes ensure all weight one 
errors, across the entire system, are detectable via error syn¬ 
dromes distinct from e*. 

DCQD involves partitioning the n-qubit system into a np- 
qubit principal subsystem P and a n a- qubit ancilla subsystem 
A. After the two subsystems are initially entangled, the princi¬ 
pal system is subjected to dynamics £ while the ancilla system 
has previously been assumed to remain isolated from any such 
process. Syndrome measurements project the state into one of 
the 2 r subspaces defined by the QEC code, and the rela¬ 
tive frequency with which each syndrome is observed charac¬ 
terizes the process matrix A schematic circuit expressing 
this partition for DCQD is shown in Fig.[lja). 


A. Clean Ancilla DCQD 

The probability for a state p = \^){^\ to be projectively 
measured into the subspace ^ with the error syndrome is 

Pi = Tr [Uip\ (3) 

where 11^ is the projector onto We assume a one¬ 

dimensional codespace and write the projector into each or¬ 
thogonal subspace in the stabilizer basis as = |i)(i| where 
| i) = Ei |0) for the unique operator Ei e E. (Our results will 
also hold when the syndrome subspaces are j = 2 k dimen¬ 
sional and the subspace projectors are Hi = JV \ij)(ij\-) Let 
the encoded state of the system be initialized as pL = |0)(0|, 


















3 


Initialize 

Noise 

Apply 

Measure 

Determine 

1 

^ J 

£ - 

Oj 

- ei- 

Xmn(Oj) 

J 


FIG. 2. Schematic of the DCQD process. Begin with the state 
p = |0)(0| initialized in the codespace of S\. Next, subject p 
to some dynamics resulting in £(p). A pre-processing operation 
Oj — {E, Uj , Pj } is then applied just prior to the stabilizer gen¬ 
erators yielding an error syndrome e%. An ensemble of syndrome 
measurements given the pre-processing Oj is used to deduce part of 
the process matrix Xmn(Oj) — {x™, Imxji, Re^Jz} as described 
in Eqs. [4l5]6| 


so that before any quantum operation is applied the trivial syn¬ 
drome eo occurs with probability po = 1. 

After £ acts on P, the probability for each error syndrome 
becomes EE) 


Pi = Tr 


n,£ 

rnn 


= Xu 


(4) 


where we use the QEC condition in Eq. ([2]) and perform the 
trace in the stabilizer basis. Therefore, the quantum dynam¬ 
ical populations (diagonal elements) of x are simply the rel¬ 
ative frequencies with which each syndrome appears in an 
ensemble of stabilizer generators measurements. 

Off-diagonals of x represent quantum dynamical coherence 
and are similarly measured by first applying the unitary oper¬ 
ation Uj = (1 + iFj)/y/2 , where Fj E E is a member of 
the Pauli group with a trivial Pauli phase factor of +1. The 
corresponding probability for each syndrome measurement is 
then 


Pi{Uj) = Tr [n iUj£{p)U] 
Xu + xjj 


(5) 


- Im (4>jXJi) 


where </>j E {±1, =bi} is a Pauli phase factor which, along 
with Fj, depends on the indices i, j according to (j)jFj = 

fUj m. 

Applying the unitary Uj enables measurements to probe ei¬ 
ther the real or imaginary part of Xij • The complimentary part 
of a given element Xij is recovered by applying instead the 
projective measurement Pj± = 1 dfc Fj. This corresponds to 
measuring an eigenvalue of ±1 for the operator Fj. In this 
case the syndrome probabilities become 


Pi(Pj±) = Tr UiP j± £(p)P] ± 
Xu + xjj 


± Re ( j>jXJi ) 


( 6 ) 


for the ±1 eigenvalue. 

Equations (4)-(5) represent a system of linear equations that 
determine the elements Xij °f the process £. Direct charac¬ 
terization can either be used to construct the complete process 
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TABLE I. Error syndromes e* for states \i) — Ei\ 0) indexed by 
the integer i for the group of located errors Ei E Ep . As evident 
through the one-to-one correspondence between the located error op¬ 
erators (Ei) and the syndromes (a), the code Si is non-degenerate 
with respect to located errors Errors involving weight-one ancilla op¬ 
erators (Ep 0 Ea E E) are associated with error syndromes whose 
first two values are either 01,10 or 11. These errors are filtered out 
and do not affect the constructed Xmn- 


matrix or it may be applied partially to characterize only spe¬ 
cific elements of x- Partial characterization is especially use¬ 
ful when a priori knowledge about the quantum dynamics is 
available. For example, partial dynamics can determine the 
relaxation times Ti, T 2 efficiently by only characterizing xtz 
and xxv 0. 


B. [[6,0,2]] Concatenated Ancilla Code and Error Filtering 

In its current form, DCQD assumes the ancilla system is 
noiseless and therefore any non-trivial syndrome measure¬ 
ments are attributed to the process acting solely on the prin¬ 
cipal system. This provides a justification for interpreting the 
syndrome statistics in terms of the quantum process defined 
by Eq. 0 acting on the principal system. Schemes assuming 
noiseless A require the operator set E have support only on P 
in order to decode the syndrome. (§]1TZ ] Relaxing the assump¬ 
tion of ideal ancilla would introduce ambiguity into syndrome 
interpretation. For example, in the Bell state previously used 
for DCQD E, a stabilizer measurement cannot discriminate 
between a process that invokes no error on A and a nontriv¬ 
ial error on P and one that induces an error on A while P is 
unaffected. 

In practice, this ambiguity leads to errors in the characteri¬ 
zation of the principal system as realistic ancilla also undergo 
quantum process. Note that the error ambiguity, seen for ex¬ 
ample in Bell state DCDQ 0, comes from the invariance of 
the error syndrome under the interchange A o P. To re¬ 
solve this ambiguity, and differentiate between the different 
physical scenarios that lead to the same syndrome, we con¬ 
catenate the ancilla A qubits. Concatenation invalidates the 
mapping A o P, thus removing the syndrome ambiguity. 
Concatenated ancilla qubits involve additional stabilizer gen¬ 
erators such that the code detects low weight processes ex¬ 
clusive to system A. These new syndromes can be used to 
filter the characterization data by rejecting those values that 
indicate errors on the ancilla. This offers an improvement to 
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a fundamental limitation of code-based process tomography. 
Moreover, the concatenated ancilla further partition the set E 
into a set of located errors with support on only the princi¬ 
pal system and a set of unlocated errors whose support is the 
composite system. 

We now outline our main result, the construction of a code 
characterizing a np = 2 qubit principal subsystem with noisy 
ancilla. The characterization code must first satisfy the located 
quantum Hamming bound J2j =o (™ P )^ k < 2 n l20l . The k = 
0 located Hamming bound (we choose k = 0 since we wish to 
minimize overhead and are note interested in encoding logical 
information) is saturated for n = 4, so we use the [[4, 0, 2]] 
code So = (XIXI, IXIX, ZIZI, IZIZ) to characterize P. 
However, as discussed in the last section, this code leads to 
the mistaken interpretation that processes acting on A (qubits 
3,4) characterize P (qubits 1,2) under the interchange A <G> P 
as is obvious from the symmetry of the generators. 

We can remove the syndrome degeneracy by encoding each 
physical qubit in A with a second error detection code. This 
form of code concatenation enables the detection of processes 
that occur on only the ancilla. Concatenation of the ancilla is 
also compatible with DCQD, as the first level stabilizers So 
are still used for direct process characterization. A schematic 
of this process is shown in Fig. [T] The additional resources 
required for encoding the ancilla can be managed by adjusting 
the error detection properties of the second code. We encode 
the two ancilla qubits from the original characterization code 
with a [[4, 2, 2]] code that is capable of detecting weight-one 
operators G3ED. This brings the total number of qubits to 
six and forms a [[6, 0, 2]] code. 

The encoding [[4,2,2]] stabilizer group is 5^ = 

{XXXX,ZZZZ), where we choose X x = XXII, Z 1 = 
ZIZI,X 2 = IXIX,Z 2 = IIZZ as representative logi¬ 
cal operators for the two encoded qubits. Ancilla concate¬ 
nation means replacing the ancilla qubits in So with the logi¬ 
cal qubits from Se , that is, X(Z) 3?4 X(Z) i j2 and So ^ 
{XIXI, IXIX, ZIZI, IZIZ). Expressed in terms of its 
generators, the newly formed code is 

Si = {IIXXXX, IIZZZZ, XIXXII, (7) 

ZIZIZI, IXIXIX, IZIIZZ). 

It is clear from Eq. 0 that either one or both of the first 
two generators anti-commute with all weight-one errors on 
the four qubit ancilla subsystem. The remaining generators 
associate a unique error syndrome to all located errors on P 
when A is noisless. Additionally, as detailed below, this code 
detects all errors occurring on the first two qubits P simulta¬ 
neous to any weight-one errors on A. 

The group Si stabilizes the (unnormalized) one dimen¬ 
sional codespace 

| 0 ) = 1000000 ) + 1001111 ) + 1010101 ) + 1011010 ) ( 8 ) 
+ 1100011 ) + 1101100 ) + 1110110 ) + 1111001 ). 

In addition, S\ partitions the Hilbert space into 64 one¬ 
dimensional orthonormal subspaces Ml, each of which is 
identified by a unique error syndrome e t . 


We now detail the set of errors E for which QEC condition 
in Eq. ^ is satisfied. We begin by structuring E into two dis¬ 
joint sets based on location of the induced errors or process. 
The first set consists of located errors acting on the principal 
system P, i.e., operators of the form cr 1 cr 3 Till. We denote 
this set of 16 errors, which forms the Pauli two-qubit group 
modulus phases, by Ep = V 2 /{±i,±l}. It shall also be 
useful to refer to the 12 possible weight-one ancilla errors as 
Ea = {Xi,Yi, Zi} where i = 3,4, 5, 6 are the sites compris¬ 
ing A. The second set of errors, with 192 elements, consists 
of the tensor product of located and ancilla errors Ep (8) Ea . 
Using the above definitions, the set of detectable processes is 
the disjoint union E = Ep 0 1a U Ep 0 Ea- 

The noisy ancilla filtering properties of this code are evident 
upon inspecting the syndromes pertaining to the sets Ep 01 a 
and Ep 0 Ea . The code is non-degenerate for the set of lo¬ 
cated processes, i.e., choosing E a ,Ei) G Ep the code matrix 
in Eq. ^ becomes C a b = S^ for the state in Eq. ([5]). El¬ 
ements Ei G Ep map the codestate to distinct orthogonal 
states | i) = £^|0), with distinct syndromes for i G [0,15]. 
This group of located processes commutes with the first two 
generators in Eq. 0- so the corresponding syndromes are of 
the form = (0,0, e* 3 , e^, e^, e^). Table [I] enumerates the 
syndromes associated with all located errors. Syndromes that 
begin with “00” indicate A is error free and that the corre¬ 
sponding measurement is accurate for characterizing x as de¬ 
scribed in Eqs. |4|5|6~ 

The code is degenerate for processes E a , E^ G Ep 0 Ea- 
This result is expected since there are 2 6 = 64 syndromes 
and 208 operator elements in the set E. The remaining 192 
processes Ej G Ep 0 Ea (j G [16, 207]) map the codeword 
onto the remaining 48 orthogonal states, | i) with i G [16, 63]. 
Those processes that have an odd weight support on A anti¬ 
commute with either one or both of the first two generators 
of Si. Consequently, syndromes that begin with 01,10,11 
indicate that noise was detected on the ancilla. Because these 
syndromes are degenerate we cannot know exactly which pro¬ 
cess corrupted the ancilla. Therefore, this data is filtered out 
from characterizing the principal system. 


III. NUMERICAL CHARACTERIZATION OF 
AMPLITUDE DAMPING CHANNEL 

We now test the procedure outlined in the previous sec¬ 
tion (and Fig. [2]) by numerically constructing the process ma¬ 
trix for the well known amplitude damping (AD) channel. 
While we could characterize arbitrary noise on P, for clar¬ 
ity we consider the case when only the first qubit experi¬ 
ences AD. The principal system process matrix is £p D (p) = 
^2 E^ D pE^ D ^ where AD channel is written in terms of 
the Kraus operators Eq D = (1 + y/1 — 7)1/2 + (1 — 
yT - 'y)Z 1 /2,E^ D = y/i(Xi + iYi)/2. Expressing the 
X matrix in the Pauli basis Fi = {I, X,Y, Z} the only 
non zero process matrix elements appear along the diagonals 
Xii = (1 + VI _ l) 2 /^Xxx = Xyy = 7/4,Xzz = (1 - 
VI - 7) 2 / 4 and anti-diagonals xiz = Xzi = x/^Xyx = 
~Xxy = n/ 4 (see Fig. 0- To test our code in the pres- 
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FIG. 3. Simulated AD channel process matrices constructed from 
ensembles of syndro me me asurements according to the DCQD pro 
cedure, namely Eqs. EHf 


Probabilities are determined from 10 


,6 


Monte-Carlo events using the numerical parameters 7 = .4, p = .1 
were used for the amplitude damping and depolarizing channels re¬ 
spectively. The real and imaginary parts for x constructed with a 
noiseless ancilla is given in panels (a,b) and its difference from the 
theoretical value appears in panel (c). A noisy ancilla reduces the 
accuracy of the DCQD procedure as seen in panels (d-f) for which 
a standard [[4, 0, 2]] has been used. As seen in panels (g-i), weight- 
one ancilla errors are mitigated by utilizing a concatenated [[ 6 , 0 , 2 ]] 
code. The constructed x matrix is characterized by a high degree of 
fidelity, F{££ D (p),£ti 6 ’°’ 2 U(p)) = .9884. 


on qubit 1 with a noisy A implemented with detection being 
done by a non-concatenated [[4,0,2]] DCQD code and (iii) 
AD on qubit 1 with a noisy A using the [[6, 0, 2]] code given 
in Eq.[7]to determine x* 

The process matrix constructed in scenario (i), i.e. for a 
noiseless ancilla system, is shown in Fig.[3]panels (a,b) and is 
compared to the theoretical result in panel (c). Finite sampling 
causes a small discrepancy between the theoretical and the 
simulated result as seen in panel (c). Next, we simulate case 
(ii) involving the four qubit non-concatenated DCQD code in 
which every possible syndrome is used to determine the x* 
The absence of a filtering process means that each error occur¬ 
ring on the ancilla system corrupts the syndrome probabilities 
which, in turn, determine Xi,j- The simulated x matrix is pre¬ 
sented in Fig.[3]panels (d,e) and its distance from the clean x is 
given in panel (f). Finally, for case (iii), we construct x using 
the ancilla concatenated code (Eq|7} and present the results in 
panels (g-i). Inspecting panels (c), (f), and (i) it is obvious that 
the x matrix constructed with the concatenated code is more 
accurate than the standard non-concatenated code. 

To quantify this difference we calculate the fidelity,defined 
as F(p,cr) = Tr [^/\/po^/p\ for two density matrices 
p, a, between the theoretical and numerically constructed 
X coefficients. The states we calculate the fidelity of 
are one qubit states subjected to our constructed ampli¬ 
tude damped channels and the theoretical channel, that 

is: «***(,) = E«•***« = 
Em„ Xrnn ,2]] F rn pFl, and £p D (p). Using an initial single 
qubit state p = |0)(0| we find F(£p D (p),£^ 6 ’ 0,2 ^(p)) = 
.9884 and E(£p D (p), £ it 4 ’ 0 ’ 2 !] (p)) = .9165 which represents 
a 10% improvement in the fidelity for the constructed process 
matrix for the specific case of p = 0.1. 


ence of a noisy A subsystem we take the state £p D (p) and 
subject it to an additional depolarizing (DP) channel acting 
independently on each ancilla qubit. We construct the chan¬ 
nel via a composition of single qubit DP channels so that 

£ a P (p) = £ 3 P (p) ° £ 4 P (p) ° £ 5 P (p) ° £ 6 P (p) where 
£?■ p (p) = (1 - p)p + p{X iP Xi + Y iP Yi + Z iP Zi)/ 3 and 
/(p) o p(p) = f(g(p)) denotes the usual functional composi¬ 
tion of mappings. The resulting state is £^ p (£p D (p)) where 
for order of the the independent P, A channels is arbitrary. 

To simulate experimental measurement statistics we per¬ 
form a Monte-Carlo simulation in which we project the state 
£j± p (£p D {p)) into the ±1 eigenstate of each generator in 
Eq. [7] with probability Tr [(1 ± gi)£^ p (£p D (p))]. This pro¬ 
cedure is repeated for all six syndromes with the ±1 eigenval¬ 
ues for each generator defining a single measured syndrome 
generator. The probabilities for each clean syndrome, i.e., 
“00” syndromes, with respect to all clean results is used to 
determine the x elements by Eqs. 0, 0 , and ([6]). Following 
this procedure we perform a Monte-Carlo simulation of the 
following three scenarios: (i) the AD channel £p D (p) acting 
on qubit 1 with a noiseless ancilla A (£j^ p (p) = p), (ii) AD 


IV. FILTERING FAILURE RATE 

Figure [3] shows how encoded ancilla can improve pro¬ 
cess characterization as compared to previous QEC-based 
schemes. This works by detecting noisy ancilla operations 
and filtering out those measurements, thus improving the ac¬ 
curacy of the constructed process matrix. We can explain the 
improvements in characterization fidelity in terms of a gain 
in the signal to noise ration for the characterization process. 
In particular, we define the signal as the constructed process 
matrix elements, which are directly related to a measured syn¬ 
drome e t . Noise is then those syndromes in which an error 
Ej occurs on P but we measure a syndrome with j ^ i. 
The rate at which this occurs is represented by pp in Fig. [4] 
To quantify the improvement in the signal to noise ratio we 
now analyze a simple model, the DP channel £^ p (p), and 
compare the probability of failure for concatenated and non- 
concatenated characterization codes. 

Using the definition of noise, we say that the filter fails if 
data collected from a noisy event is used to characterize Xmn- 
In the [[6, 0, 2]] code, failure cannot be due to any weight-one 
errors on A since they, and their product with all located errors 
(Ep (8) Ea), are within the detectable errors set E. The filter 
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FIG. 4. Probabilities of located error syndrome to occur in the pres¬ 
ence of a noiseless principal system P in the presence of independent 
depolarizing channels acting on the ancilla A with probability p. The 
blue dashed line (Vt = (1 — p) 4 ) shows the probability for the iden¬ 
tity operator to occur while the probability for the identity syndrome 
to appear in the Monte-Carlo simulation is given by the blue circles 
(denoted by pt). Stabilizer operators therefore occur with probabil¬ 
ity Apt = pt — Vt- Orange circles (poo) denote the rate at which 
the remaining located (“00”) syndromes occur. Green circles (dashed 
lines) denote the simulated (theoretical) failure rate pf — Poo + Api 
(Vt = P 2/3 + 2p 3 / 9 + 21p 4 /81, see text for derivation) indicating a 
located error syndrome from Tab.[I]is caused by an operator different 
than Ei . 


does however fail in the presence of some weight-two errors 
which commute with the first two generators in Eq. [7] In the 
DP channel f£ p (p) the probability for a weight j to occur is 
Pj = (1 — p) 4- - 7 (f)pi • Notably, the probability for the weight 
0 “error” 1 to occur is the probability that the identity occur 
on each qubit po = (1 — p) 4 which appears as the dashed blue 
line labeled V\ in Fig. [4] Ancilla errors outside the correctable 
error set Ea occur with probability p> 2 = P2 +p 3 +P4. How¬ 
ever, not all of the errors with weight > 2 will lead to faulty 
characterization data since many most of them will still lead to 
syndromes beginning with one of 01,10,11 and therefore do 
not corrupt the constructed In these cases we discard the 
data point because it is (correctly) assumed that some error 
has occurred on A. 

To confirm our estimates, we numerically calculate the fail¬ 
ure rate with 10 6 Monte-Carlo simulations of the composite 
depolarizing channel f£ p (p). With a single exception, the 
failure probability is by definition the number of syndromes 
beginning with 00 divided by the total number of randomly 
generated errors. The exception comes from the ambiguity of 
whether the syndrome eo = {0, 0, 0, 0, 0,0} should count to¬ 
wards the error rate, as eo may be generated by the identity 
mapping or by any element in the normalizer group 
i.e. the group of errors commuting the all stabilizer elements. 

However, we know that the identity operator (l 04 ) occurs 
with probability V\ = (1 — p) 4 as illustrated by the blue 
dashed curve in Fig. [4] We determine the rate for erroneous 
identity-like syndromes to be A p t = Pt — T\, the differ¬ 
ence between V\ and the numerical rate at which we mea¬ 
sure the identity syndrome (blue circles in Fig. [4]). In Fig. [4] 
the green circles represent the total failure rate obtained by 


adding the identity probabilities difference to the probabil¬ 
ity with which all other located syndromes occur. Enumer¬ 
ating the number of weight 2,3, and 4 errors which com¬ 
mute with the first two generators of S\ and the probability 
with which they occur we find the probability of failure to be 
Vf = 2 p 2/3 + 2p 3 /9 + 21p 4 /81 where P 2 , 3,4 is the for prob¬ 
ability for an error of weight 2,3, or 4 to occur. This function 
of Vf is plotted as the dashed green line in Fig.[4]and exactly 
matches our numerical data. The leading term in Vf goes 
as 0(p 2 ) in contrast to to non-concatenated DCQD schemes 
whose failure rate gores as O(p), the probability for weight- 
one errors, explains the sharp contrast in the constructed pro¬ 
cess matrices in the second and third rows of Fig. [3] 


V. DISCUSSION 

We have introduced a DCQD code that directly character¬ 
izes the quantum dynamics of a principal system with as¬ 
sistance from a noisy ancilla system. Within the stabilizer 
framework, we show that ancilla noise can be distinguished 
from processes acting on the principal system by using syn¬ 
drome value as a filter for non-trivial ancilla processes. For 
the example of DCQD with a [[4,2,2,]], we have concate¬ 
nated the ancilla qubits for purposes of detecting weight-one 
processes, and compared the characterization of an amplitude 
damping process on the principal system using three different 
approaches: (i) clean ancilla system, (ii) noisy ancilla using 
a standard DCQD, and (iii) noisy ancilla using our concate¬ 
nated [[6,0,2]] code. Our numerical simulations found that 
the process matrix constructed using the six-qubit code shows 
a marked improvement in fidelity over the non-concatenated 
approaches. 

Our motivation for encoding the ancilla qubits has been to 
filter out those measurements that correspond to unwanted 
data. From this perspective, ancilla encoding represents a 
form of filtering the dynamics to isolate non-trivial processes 
acting only on the principal system. We have argued that fil¬ 
tering increases the signal-to-noise ratio for process character¬ 
ization, as measured by the gain in fidelity of the constructed 
matrix. Of course, the gain for process characterization de¬ 
pends strongly on the details of the ancilla filter. For exam¬ 
ple, the 6 -qubit code introduced here detects only weight-one 
ancilla errors and their product with located principal system 
errors Ep ( 8 ) Ea- When higher weight errors are common, 
the benefit of this ancilla encoding diminishes, and larger dis¬ 
tance codes are needed to filter higher weights processes. For 
example, a distance 4 code that detects all weight-2 ancilla er¬ 
rors will have a filter failure rate that scale as 0(p 3 ) with p the 
ancilla error rate. We could also have used a non-degenerate 
[[5,1,3]] code to encode the ancilla, where each detectable 
error would have a unique error syndrome. In this case, each 
syndrome would be used without a filtering procedure. In gen¬ 
eral, one can improve the signal to noise ratio at the expense 
of additional ancilla qubits and larger codes. 

Additionally, we have taken k = 0 throughout thought 
this work, but we could have used a k / 0 code satisfying 
a generalized Hamming bound l20l . For example, a non- 
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concatenated code performing error correction on two qubits 
with another encoded is provided in Ref.lT2l Equations ([4])-([6]) 
are easily generalized using the higher dimensional projectors 
Pii resulting in a code which detects ancilla errors while en¬ 
coding some non-trivial quantum information. 

Recent progress in realizing stabilizer QEC circuits with 9 
and 4 qubits on different lattice configurations suggest that the 
implementation of these ideas should be experimentally fea¬ 
sible in the near future [17, 18]. In particular, it is worth not¬ 
ing that the characterization processes described here and in 
earlier DCQD work do not require active, feed-forward error 
correction for purposes of implementation. Consequently, the 
use of QEC-based DCQD appears to be a natural way point 
toward the demonstration of error-corrected computation. 
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